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Abstract 

The notion of the optimized perturbation, which has been successfully ap- 
plied to energy eigenvalues, is generalized to treat wave functions of quantum 
systems. The key ingredient is to construct an envelope of a set of perturba- 
tive wave functions. This leads to a condition similar to that obtained from 
the principle of minimal sensitivity. Applications of the method to quan- 
tum anharmonic oscillator and the double well potential show that uniformly 
valid wave functions with correct asymptotic behavior are obtained in the 

first-order optimized perturbation even for strong couplings. 
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Naive perturbative expansions are usually divergent asymptotic series [[!]. Thus various 
summation techniques such as the Borel and Pade methods have been developed [0]. In 
recent years, a new summation method "optimized perturbation theory (OPT)" has been 
studied. This combines straightforward perturbation with a variational principle and im- 
proves even the Borel non-summable series. Variants of OPT are called in various names, 
the (optimized) ^-expansion, the variational perturbation, the renormalized strong coupling 
expansion, and so on ||. One of the key ingredients in the method is the principle of minimal 
sensitivity (PMS) (§], which says that physical quantities calculated in perturbation theory 
should not depend on any parameters absent in the original Hamiltonian. 

The method has been applied to improve the energy eigenvalues of the quantum anhar- 
monic oscillator (AHO) and the double- well potential (DWP); AHO and DWP are the useful 
laboratories to test a new non-perturbative method. A rigorous proof has been also given for 
the convergence properties of the perturbation series for energy eigenvalues of AHO/DWP 
0. The method is also being extensively applied to perturbation series for energy and 
partition functions in quantum field theories |J. 

To the best of our knowledge, the wave functions has never been explored in OPT even for 
simple quantum mechanical problems. Needless to say, wave functions have all information 
of the system in question, including the energy eigenvalues. The purpose of this Letter is to 
show a novel generalization of the idea of the optimized perturbation for studying the wave 
functions. 

Let's take the following Schrodinger equation for one dimensional AHO/DWP; 

\{-^ ± x 2 + Ax 4 )^ n (x) = E n (\)ijj n (x). (1) 

Naive Rayleigh-Schrodinger (RS) expansion for eigenvalues E(\) = J2 n A n C n gives an asymp- 
totic series with coefficients growing as C n ~ n\ A scaled Hamiltonian in the simplest 
version of OPT (sometimes called the linear (^-expansion) reads H§ = Hq(Q) + 5 ■ Hi(Q, A) 
with 

JMn) = l(-^ + fiV), (2) 
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H^n, A) = i ((±1 - ft 2 )x 2 + Ax 4 ) . (3) 

Hs interpolate the free Hamiltonian with a frequency Q and the full Hamiltonian; Hs=o = 
Hq(Q) and Hg=i = H. The essential idea of OPT is to choose a suitable Q so that the 
RS perturbation E(X) = J2n^ n C n becomes a convergent series with errors exponentially 
suppressed for large n. PMS in this case is one of the possible criteria for choosing Q, 
although it is not a unique one [[| . 

To show how the naive application of the idea of OPT to the wave function fails, we first 
expand the wave function of H$ by the eigenfunctions of H , 

i> k &ty = J2c jk (n)rf(x-,n). (4) 
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Cjk(£l) is evaluated in the RS perturbation theory with 5 being an expansion parameter, 

CjkiP) = Cj k + 5 ■ Cj k H . With a suitable Q, one might expect ip(x, ft) at 8 = 1 converges 

rapidly to the exact wave function ip exact (x) as in the case of the energy eigenvalues E. 
However, this is not the case as long as Q is independent of x. In fact, in any finite orders 
of the expansion in (f|), ip n (x — ► oo; Q) ~ exp(— Qx 2 /2), while the exact wave function in 
(H) behaves as 

^(ar-ooJ^ejqpC-^llarl 3 ). (5) 

Thus, the series (|j) can never reproduce the true asymptotic behavior no matter what Q is 
chosen as long as Q is constant. 

To find a way out, let us describe the situation in geometrical terms. The perturbative 
wave function (^) with constant Q gives a locally valid approximation around certain x 
where the perturbation 5 ■ Hi(Q) is sufficiently small; however, this is not a uniformly valid 
solution in the entire domain of x. Conversely speaking, if Q is varied, we have a family 
of functions {ip(x; parametrized with Q, each of which is locally valid around certain 
x. Therefore a uniformly valid wave function may be constructed as an envelope of the 
family of perturbative ones. To obtain an equation to determine the envelope, we recall the 
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classical theory of envelopes ||: Suppose one has a set of functions f(x;a) parametrized 
by a. Then the envelope of {f(x;a)} a can be obtained by solving d f (x; a) / da = and 
substituting the solution a(x) into the original function i.e. f(x; a(x)). In our case, f = if) 
and a = Q. Thus we reach a condition, 

dip(x; f2) 



dtt 



= 0. (6) 

5=1 

This condition has another interpretation in term of PMS: The exact wave functions of (|]) 
do not depend on the fictitious PMS parameter Q at 5 — 1. So one would impose the 
condition that even the approximate wave function should not depend on Q, which leads to 
eq.(|6]). A similarity of this idea to that of the renormalization group a la Gell-Mann-Low is 
also to be noted. 

At this point, one may ask whether ip(x;Q(x)) thus obtained satisfies the original 
Schrodinger equation where d/dx acts on Q(x) too: (d/dx)ip(x; Q(x)) = (d/dx)ip(x; fl(x)) + 
(dQ(x) / dx) ■ (8/ dfl)if)(x; Q(x)). The answer is yes, since the second term vanishes on account 
of the condition (||). The same is true for d 2 /dx 2 . Hence the improved wave function satisfies 
the Schrodinger equation at every point x at least up to the same order of the (^-expansion 
employed. 

Solution of (|) gives Q as a function of x (and A). For AHO/DWP, eq. (|6|) turns out to 
be an algebraic one. The envelope wave function ip(x;Q(x)) obtained by substituting the 
solution of into (|j) turns out to be a uniformly valid approximation for ^ cxact (x) in a 
global domain even in the lowest order of the ^-expansion. 

Let us demonstrate the above points explicitly for AHO/DWP with the first order 5- 
expansion in which Cjk = Sjk — (j \ Hj \ k)/(E® — E®). The condition eq. (H) gives an 
algebraic equation for fl as 

fi(±l - Q 2 ) (A n H n (VQx) + An{n - l)A n . 2 H n ^ 2 (Vnx)) (7) 
+A • (B n H n (VUx) + 4n(n - l)B n ^ 2 H n ^ 2 (Vnx)) = 0, 



where H n (z) is the n-th Hermite polynomial, and the functions A n (y/Hx) and B n (\/VLx) 



are given as follows; A n (z) = lQz 4 - 16(2n - l)z 2 - A(8n + 7), A n _ 2 (z) = -16, B n (z) = 
8z 6 + 4(2n + 9)z 4 -2(24n 2 + 4n-33> 2 -9(8n 2 + 26n + ll), and B n _ 2 (z) = -2(8z 2 + 18n + 9). 
Several remarks on Eq.(^) are in order here. 

1. Eq.([?D gives an algebraic equation of degree m + 5 for Q with n = 2m or 2m + 1. As in the 
usual PMS condition applied to the energy eigenvalues, one needs some principle to select 
a relevant root. The true solution Q(x) is selected so that the local energy £ n (x) defined by 
ip(x; fl(x))^ 1 H if)(x; fl(x)) is (roughly) constant as a function of x. 

2. For A = 0, there is an obvious solution Q(x) = 1 for (|7|). Other solutions do not satisfies 
S n (x) ~constant. 

3. For x — > oo, eq.(0) with an ansatz VL{x) ~ Cx p (p > 0) gives an equation independent of 
n. It has a solution C 2 = A/2 and p = 1. Thus the asymptotic form reads 

^(x;fi(x))^exp(-^/||a;| 3 ) (8) 

which has a same asymptotic behavior with (0). This feature is never obtained by any 
finite-order RS perturbation as we mentioned. A slightly different numerical coefficient of 
| x\ 3 is to be improved with higher order ^-expansion 0. 

4. Although ip(x; Q) with constant Q is normalized as Pf^, | ip(x; Q) \ 2 dx = 1, the norm of 
i/j(x;Q(x)) deviates from unity due to the x dependence of Q. In the following figures, we 
will always show renormalized wave functions. 

In fig. |l|, the real and positive solutions of (|7|) for the ground state of AHO with inter- 
mediate coupling (A = 1) are shown as a function of x. The local energy is approximately 
constant for the branch I. In the interval 0.684 < x < 0.780, Q(x) becomes complex with a 
small imaginary part. When we construct a continuous wave function in this interval, we 
adopt a real part of Q(x) as a simplest way of interpolation fllTf . We expect that the width 
of this interval and the associated small imaginary part of Q(x) will vanish in the higher 
order 5-expansion. This must be, however, checked explicitly, which is now under progress. 

In fig. 0, the ratio R = ip(x; VL(x)) / if) eyiact (x) as well as ip(x; Q(x)) for the ground state of 
AHO are shown for intermediate coupling (A = 1) and strong coupling (A = 100). ip exa " ct (x) 



is numerically solved in double precision using the shooting method with FORTRAN [12| . 

It is not until x exceeds x c = 2.725(1.211) for A = 1(100) that the relative error of 
ijj(x; Q(x)) to the exact one becomes more than 10%. Such a value of x c is large enough in 
the sense that the absolute value of the wave function at x c is already as tiny as ip(x c ) = 

0. 00016.0.0012). Hence our first order results can be used as an exact one for practical use 
of the wave function. In fig. [|, 0(1) deviation is seen for quite large values of x. This is 
due to the approximate asymptotic behavior of ip(x; Q(x)). The ratio R for large x can be 
estimated from (||) and (f|) as R ~ exp (-omVx I X . 

To see the accuracy of the wave function more quantitatively, a comparison of ip exa,ct (x) 
and i/j(x; Q(x)) for A = 100 is shown in Table | for the ground state of AHO. 

As we have mentioned before, ip(x; Q(x)) is the envelop of the perturbative wave functions 
parametrized by Q. In fig. |3], thin lines are the wave functions for A = 1 with constant Q 
varied. There are two envelopes of this family of wave functions, which are labeled by I and 
II in the figure. ip(x; Q(x)) is obtained by the envelop I. 

In Table [Tj]> expectation values of several operators including the full Hamiltonian H with 
respect to ip(x; Q(x)) are compared with the exact numbers. The agreement is excellent as 
a first-order approximation. 

Some final remarks are in order here. 

1. Although we have shown only the ground state wave function because of the limitation 
of space, we have also studied the excited states {n ^ 0) and found that i]){x] &{x)) are as 
good approximation for ^ exact ( 

the n = case. 

2. The present method here can be used to study the double well potential where the sign of 
the x 2 term in ([]]) is negative. In the first order 5-expansion, we have found that the exact 
wave function is accurately reproduced for strong couplings, but the approximation breaks 
down qualitatively for A < 0.065. In such weak couplings, or equivalently a high barrier at 
the origin, the ground state wave function is well localized in the two degenerate minima 
of the double well potential. Therefore, it is natural to expect that the naive (^-expansion 
scheme using Hq(Q) with a single harmonic potential does not work. An introduction of a 
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new PMS parameter is under investigation to widen the applicability of the method. 

3. We have also performed the second order 5-expansion for AHO. The formula is more 
lengthy and the degree of the algebraic equation is higher. The detailed comparison with 
the first order result is under progress. 

4. Applications to other quantum mechanical system such as the scattering problem, Zeeman 
effect, and Stark effect [jlj are interesting. Also, applications to field theory such as the 
improvement of the space-time (or momentum) dependence of Green's functions are closely 
related to the problem studied in this work. 

5. The notion of envelopes has been found useful for global asymptotic analysis of some of 



ordinary/partial differential equations having secular terms in perturbative solutions |T0 



The present work together with those in |T0[ shows that the notion of envelopes can be a 
useful guide for the improvement of perturbative expansions and asymptotic analyses. |13 



In summary, we have developed an optimized perturbation theory for wave functions for 
AHO and DWP. The optimized wave functions are constructed as an envelope of perturbative 
ones, and the principle of minimal sensitivity can be regarded as a condition for constructing 
envelopes. The optimized wave function thus obtained in the first-order 5-expansion has a 
correct asymptotic behavior at x — > oo and is uniformly valid for entire domain of x even 
for the strong coupling. 
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FIGURES 

FIG. 1. Q(x) as roots of the algebraic equation for the ground state with intermediate coupling 
A = 1. The branch I gives approximately constant local-energy. 

FIG. 2. The ratio R = ip(x; £l(x))/ip exact (x) as well as ip(x;£l(x)) itself for the ground state. 
The solid (dashed) line is for the intermediate (strong) coupling A = 1(100). 

FIG. 3. Family of perturbative wave functions with constant Q with 0.1 < < 100 for the 
ground state with A = 1. Thick lines labeled as I and II are their envelops. 
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TABLES 





x = 0.0 


x = 0.5 


x = 1.0 


x = 1.5 


x = 2.0 


ip(x; Q(x)) 


1.170E0 


5.519E-1 


1.885E-2 


3.925E-6 


3.095E-13 


ijj exact (x) 


1.167E0 


5.564E-1 


1.979E-2 


5.168E-6 


7.982E-13 



TABLE I. Ground state wave functions of AHO for A = 100; ip(x; Q(x)) is an approximate 
wave function in the first order 5-expansion. 





(P 2 ) 


(x 2 ) 


(x 4 ) 


(H) 


A = 1 PMS 


0.83146 


0.30391 


0.25703 


0.69620 


Exact 


0.82630 


0.30581 


0.26024 


0.69618 


A = 100 PMS 


3.34430 


0.07641 


0.01579 


2.50003 


Exact 


3.30717 


0.07731 


0.01615 


2.49971 



TABLE II. Expectation values of p 2 = (—id/dx) 2 , x 2 , x 4 and the full Hamiltonian H for the 
ground state in AHO. 
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Fig.2 




